Local density of states subject to finite impurity concentration in graphene 
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It is demonstrated that there is a characteristic impurity concentration, at which variation with 
concentration and overall appearance of the local density of states at the impurity site in graphene 
are changing their behavior. Features that are prominent in the local density of states for the single 
impurity are disappearing from it when impurity concentration far exceeds this critical value. The 
impurity subsystem not only induces the rearrangement of the electron spectrum in graphene, but 
also undergoes a substantial spectral transformation by itself, which can be observed experimentally. 
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Introduction. The spectrum rearrangement, which is 
inherent in disordered systems, have been studied for a 
large variety of systems over a period of last decades 
IH El In essence, this phenomenon is grounded on the 
indirect interaction of impurities through the host sys- 
tem. The impurity state that is formed close to the band 
edge, or in the vicinity of any other van Hove singularity, 
has an effective radius that is large compared to the lat- 
tice constant. Thus, the spatial overlap of such impurity 
states occurs at the low impurity concentration (c <C 1). 
This overlap results in a radical alteration of spectral 
properties of the disordered system, which is especially 
pronounced near the energy of the impurity state. 

Predominantly, the spectrum rearrangement is talked 
about with regard to a remarkable change in the disper- 
sion relation or in the total density of states. However, 
with recent advances in the scanning tunneling spec- 
troscopy (STS), it became possible to observe directly 
the local density of states (LDOS) 0, S @|. Thus, it 
seems reasonable to inquire how the spectrum rearrange- 
ment might be reflected in the shape of the LDOS. 

Below we would try to examine whether an increase 
in the impurity concentration influence the behavior of 
the LDOS to such an extent, so that it is permissible 
to speak about the impurity induced rearrangement of 
the LDOS on its own account. It is well known that an 
increase in the impurity concentration throws a veil over 
the LDOS features that correspond to a single impurity 
in the system 0, S @ • I n this connection, it should be 
stressed that we will be well within our rights introducing 
the concept of the impurity induced rearrangement of 
the LDOS only in a case, when qualitatively different 
regimes in the LDOS behavior can be distinguished, and 
the critical concentration that indicates the transition 
from one regime to another is closely related to the spatial 
overlap of the impurity states. 

Our treatment of this issue is based on a straightfor- 
ward model for graphene ^(J- Because graphene is a 
purely two-dimensional (2D) object, the STS measure- 
ments of the electronic LDOS in this material are de- 
cidedly natural. Moreover, possibilities of the resonance 



state formation and respective features in the single im- 
purity LDOS at the impurity site and in its neighbor- 
hood have been discussed elsewhere 0, EH E3 ■ Oppor- 
tunities to check up by future STS experiments on these 
characteristics in the near field of the impurity have been 
already addressed Q , and first reports on the STS mea- 
surements in graphene have been published recently |l5| . 

Model. For the sake of simplicity we consider a sub- 
stitutional binary alloy with a diagonal disorder in the 
tight -binding approximation. This uncomplicated model 
of a disordered system is attributed by tradition to Lif- 
shitz 0. It features absolute randomness in space dis- 
tribution of impurities. We choose the asymmetric defini- 
tion of impurity perturbation. Consequently, the on-site 
potentials are Vl with the probability c, or otherwise. 
The corresponding Hamiltonian reads, 
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where n refers to lattice cells, a enumerates sublattices, 
and c na are electron creation and annihilation op- 
erators. The summation in Eq. Q is restricted to those 
sites that are occupied by impurities. Since the impurity 
perturbation is local in its character, the main physics 
of graphene in the linear dispersion domain can be cap- 
tured by a model 2D system with a single Dirac cone in 
the spectrum. Thus, the host Hamiltonian Hq can be 
written in the following way, 

Ho = X)[/(J0^(feta(fc)+/*(*)4(feta(*)]> ( 2 ) 
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where the magnitude of the hopping parameter is chosen 
so that the bandwidth is unity, when the Brillouin zone is 
approximated with a circle. Then, the diagonal element 
of the host Green's function (GF) g = (e — iio) -1 m the 
vicinity of the Dirac point is given by |13l Il7| 

g nana {e) = 9o(e) « 2eln|e| - in |e| , |e| < 1. (3) 
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Single impurity problem. In a case, when only a single 
impurity is present in the lattice, the diagonal element of 
the GF Q — (e — H) -1 at the impurity site becomes 



g (e) = g (e)/[l - V L g (e)}, 



(4) 



For a sufficiently large impurity perturbation Vl , the re- 
spective single impurity local density of states (LDOS) 
Pimp(e) = — 7T _1 Im Go(e) manifests a resonance peak fl^j . 
This peak is analogous in its nature to the resonances de- 
scribed in unconventional superconductors |J| . Its energy 
is provided by the Lifshitz equation, 



1 = VLReg (e r ) ~ 2VLe r ln|e r | 



(5) 



The solution of this equation, i.e., the resonance energy 
e r is located above the Dirac point in the spectrum for 
perturbation Vl < 0, and vice versa. This property of the 
resonance state holds valid for any two symmetric bands 
touching each other at a certain energy. By grouping 
Eqs. © and together, obtain, 



Pim P (e) = \e\/[{l-2V L e\n\e\) 2 + (itV L ef 
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Power series expansion of the denominator in Eq. © 
about the resonance energy yields, 

Pim P (e) « \e\T 2 r /{[TrV L e r } 2 [(e- e r ) 2 +T 2 r ]}, 

T r = 7r|e r |/[2|l + ln|e r ||], (7) 

Thus, the LDOS has the Lorentz profile in the vicinity 
of the resonance, when it is located so close to the Dirac 
point, that the inequality 
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is fulfilled. As a rule, a resonance is accepted as a well- 
defined, when this condition is met. 

Finite impurity concentration. It is convenient to 
consider conditional (or weighted) GFs when dealing 
with the finite impurity concentration. The conditional 
GF with the first site occupied by an impurity, 



Vl ftimpQ j 



(9) 



can be represented through the self— energy after averag- 
ing over possible impurity distributions [la . fl9j , 
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where G = (@} is the averaged GF of the system. As 
long as the impurity concentration remains sufficiently 
small to neglect multiple occupancy corrections, the self- 
energy can be approximated by means of the well-known 
modified propagator method 



The self-energy, which is identical on both sublattices 
owing to the symmetry, is also site-diagonal within this 
approximation. In order to obtain a quantity that can be 
related to Go(e) in a single impurity problem, the diagonal 
in lattice indices element of the conditional GF should be 
properly scaled with the impurity concentration, 

c- 1 Gt PMst) * 9o(e - a)/[l V L g (e - a)}. (12) 

Thus, we get an expression that formally resembles 
Eq. (0} , in which the host GF g is replaced by the propa- 
gator of the disordered system G(e) ~ g(e — a). Finally, 
the average LDOS at the impurity site can be written 
down through the self-energy, 
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With the help of the standard substitution, 
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e — a = k exp(iy), x > 0, < <p < n, 



(13) 



(14) 



the imaginary part of the self-consistency condition in 
Eq. l(TT|) becomes, 

cVl [2\xik+ (2cp - 7r) cot <p] + 
+ [1 — Vl>c{2 lnxrcos ip — (2<p — ir) sin ip)] 2 + 
+ [V L x(2ln k sirup + (2ip - n) cos (p)} 2 = 0, (15) 

where the particular form of the host GF (see Eq. 
has been taken into account. At the given impurity per- 
turbation Vl and the given impurity concentration c, this 
equation always has two solutions for the phase <p, when 
h exceeds some threshold value. In a turn, both roots 
provide corresponding magnitudes of the energy, when 
substituted into the real part of Eq. (fill) . 

e = xcos</3 + cVl [1 — Vix(2Inxcos<^— 

— (2ip — it) sin ip)] / {[1 — Vl>c(2 In x cos <p — 

— (2ip — 7r) sine/?)] 2 + [Vi>f(21n>{:sin(/5 



(2p — 7r) cos <p)] 2 }- 



(16) 



The LDOS from Eq. (|13|) can be expressed in the same 
variables too, 



Pzoc(e) = xsixnp/(iTcVl). 



(17) 



Several examples of the LDOS at the impurity site, which 
were calculated according to Eqs. (|13|l - lfT7|l . are depicted 
in Figs, mil f° r a range of impurity concentrations. 

Rearrangement of resonance peak. In a parallel to the 
single impurity problem, the resonance energy that cor- 
responds to the given impurity concentration should be 
a solution of the equation (cp. to Eq. JSJ), 



cV l /[1- V L go(e-a)}. 



(11) 



1 = V L ReG {e r (c)). 



(18) 
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By using Eq. J3J and the substitution l|14(l , this equation 
can be recast as follows, 

1 = 2T^e r (c)[ln|e r (c)/cosip r (c)| + 

+ (7r/2 - <p r (c))tancp r (c)]. (19) 

From here on, we will assume without any loss of general- 
ity that e r > to make expressions more readable. When 
the concentration is low enough to keep f r (c) small, an 
approximate solution of Eq. (|19|1 is: 



e r (c) w e r [l +7 r <^ r (c)] 



<p r {c) "C 1. 



(20) 



Thus, because j r (see Eq. (jHJ) is the resonance width pa- 
rameter that should be small for the well-defined quasilo- 
calized state, the energy of resonance varies slowly with 
tp r (c). The concentration dependence of tp r (c) can be 
obtained from the self-consistency condition l|15|) . The 
second term in Eq. H15fl is zero by the definition of e r (c), 
the remaining two lead to the relation, 

c = -2e^(c) tanyv(c) x (21) 
x (In I e r (c) / cos ip r (c) | tan ip r (c) + (p r (c) — 7r/2) . 

As above, Eq. Q21JI considerably simplifies while the phase 
at the resonance energy ip r (c) remains small, 

c w 7re2(c)^ r (c)(l + ip r {c)hr), </>r(c) < 1. (22) 

This relation reveals a hidden physical meaning of pa- 
rameter 7 r . It clearly shows that -f r is the characteristic 
phase for the resonance state. The corresponding impu- 
rity concentration is given by Eq. I|22|) with ip r (c) = 7 r , 



1 + In e r 
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(23) 



At the critical impurity concentration c r the damping 
of the resonance induced by the disorder becomes equal 
to the damping of the single impurity state, — Imer ps 
e r (c r )(p r (c r ) ~ r r , and, respectively, the LDOS at the 
resonance energy decreases to the one half of its magni- 
tude for the single impurity problem (see Eq. 0), 



Ploc(£r(c r )) 
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Neglecting the concentration offset in the resonance 
position (|2UI) . Eq. 122|) can be easily solved for the phase, 



¥V(c) « (7r/2)( v / l + (8c/c r )- 1), 



(25) 



Consequently, the concentration dependence of the res- 
onance peak height in the LDOS at the impurity site 
immediately follows from Eqs. I|24|l and (|25|l . 



y/1 + (8C/C,.) - 1 

PiocMc)) w 7j— ^ r p lm p(e r )- 



(26) 
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FIG. 1: The LDOS at the impurity site for e r = 0.01 at 
c = An x 10" 5 , n = 1, 2 ... 7. The single impurity LDOS is 
displayed for comparison by the thick line. The peak height 
is decreasing with increasing the concentration. The inset is 
showing the bottom left hand corner at a larger scale. 
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FIG. 2: The LDOS at the impurity site for e r = 0.01 at 



(36+4n)xl0" 



1, 2 ... 7. The single impurity LDOS 



is displayed for comparison by the thick line. The peak height 
is decreasing with increasing the concentration. The inset is 
showing the bottom left hand corner at a larger scale. 



With increasing impurity concentration, this height de- 
creases linearly with concentration at the start, 

Pioc{e r (c)) w [1 - (2c/c,.)]p imp (e r ), c < c r , (27) 

but then the rate of this gradual decrease slows down, 



Pioc(£ r (c)) ~ y/c r /(2c)p im p(e r ), c>c r . (28) 

The difference between these two regimes is clearly evi- 
dent from the comparison of Fig. ^ with Fig- EI in prin- 
ciple, it should be feasible to capture the change of the 
rate at which the magnitude of the LDOS at the reso- 
nance energy varies with the impurity concentration by 
STS measurements, not to mention the halving of the 
resonance peak height. Thus, we expect that the task 
of pinpointing the critical concentration of the spectrum 
rearrangement, which is a valuable parameter of the sys- 
tem, is accessible by actual experiments. 

Rearrangement of antiresonance. Apart from the res- 
onance, there is also a noticeable dip in the LDOS shape 
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close to the Dirac point in the host system, which can be 
considered like a kind of an antiresonance. At low im- 
purity concentration the position of this dip practically 
coincides with those energy in the spectrum, at which 
ip = 7r/2. At this energy the imaginary part of the self- 
consistency condition Eq. (fTT|) reads, 



Im (j = H dlp 



2cVl>Cdip In >Cdi p 



(29) 



while the effective shift due to impurities is given by 

cV L 



£di P (c) = Re a 



(30) 



As usual, the characteristic concentration for this point 
in the spectrum should be determined by the relation 
| Reer| = | Im a \ . This condition yields, 



Cdip 



= -Aet ln|e r 



(31) 



For c <C Cdip, the magnitude of the LDOS at this specific 
energy is rapidly increasing with impurity concentration, 

Pioc{t dl p{c)) a (2A){exp[-l/(2 C y L 2 )]/(2cV L 2 )}, (32) 

and then is reaching its maximum value at c = Cdip, 

Pioc{edtp(cdt P )) = (2tt\V l \)~ 1 . (33) 

According to Eq. (|3TJ)l . the dip position is gradually dis- 
placed approximately by the amount cVl with increasing 
the impurity concentration for c <C Cdip, and, finally, this 
dip totally disappears from the LDOS at the impurity 
cite at c ~ Cdip- The described detachment of the LDOS 
curve from the energy axis can be distinctly seen in the 
inset of the Fig. [2] 

Discussion In a system with the linear dispersion, the 
GF variation with the intersite distance has the charac- 
teristic radius that is proportional to l/\e\. Since in a 2D 
system the average distance between impurities is propor- 
tional to 1/y/c, the anticipated magnitude of the critical 
concentration of the spectrum rearrangement is sitting 
around e 2 ,. Two critical concentrations ll'-'il) and (|31|l that 
were argued for above are bracketing from both sides this 
crude estimation. There is an intimate interdependence 
of this LDOS rearrangement and the rearrangement of 
the whole electron spectrum, which have been studied in 
fl3| . It is not difficult to verify that the validity criterion 
introduced in 01 i s completely supporting the results 
derived here. Within the plain impurity model adopted, 
the well-defined resonance state is possible only for a 
strong impurity perturbation. However, the restrictions 
imposed on the resonance appearance f or, say the double 
impurity should not be that severe |li|2ll|22|, while the 
main physics of the LDOS rearrangement must remain 
substantially the same. Analogous effect can be achieved 
by a tuning of the impurity-host hopping parameter. 



Conclusion. In summary, we have demonstrated that 
the concept of the LDOS rearrangement is fully justi- 
fied. The sharper and the closer to the Dirac point is the 
resonance, the lower is the critical concentration of the 
LDOS rearrangement. The respective critical concentra- 
tion should not be exceeded for the resonance peak and 
the antiresonance dip to be discernible in the shape of 
the LDOS at the impurity site. We presume that the 
LDOS rearrangement is not specific to graphene or re- 
lated systems with the Dirac dispersion, and should oc- 
cur in virtually any system that manifests impurity states 
in the close vicinity of the van Hove singularities in its 
spectrum. 
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